Loschmidt echo decay from local boundary perturbations 
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We investigate the sensitivity of the time evolution of semiclassical wave packets in two- 
dimensional chaotic billiards with respect to local perturbations of their boundaries. For this 
purpose, we address, analytically and numerically, the time decay of the Loschmidt echo (LE). 
We find the LE to decay exponentially in time, with the rate equal to the classical escape rate from 
an open billiard obtained from the original one by removing the perturbation-affected region of its 
boundary. Finally, we propose a principal scheme for the experimental observation of the LE decay. 



The study of the sensitivity of the quantum dynamics 
to perturbations of the system's Hamiltonian is one of 
the important objectives of the field of Quantum Chaos. 
An essential concept here is the Loschmidt echo (LE), 
also known as fidelity, that was first introduced by Peres 
[H and has been widely discussed in the literature since 
then [2j. The LE, M(t), is defined as an overlap of the 
quantum state e~ lHt / h \4>o) obtained from an initial state 
\4>o) in the course of its evolution through a time t under 
a Hamiltonian H, with the state e~ lHt / n \4>o) that results 
from the same initial state by evolving the latter through 
the same time, but under a perturbed Hamiltonian H 
different from H: 



M(t)= (<j>, 



D \e ' e 



iHt/hl 



(1) 



It can be also interpreted as the overlap of the initial 
state |<^o) and the state obtained by first propagating 
|(/>o) through the time t under the Hamiltonian H, and 
then through the time —t under H. The LE equals unity 
at t = 0, and typically decays further in time. 

Jalabert and Pastawski have analytically discovered Q 
that in a quantum system, with a chaotic classical coun- 
terpart, Hamiltonian perturbations (sufficiently week not 
to affect the geometry of classical trajectories, but strong 
enough to significantly modify their actions) result in the 
exponential decay of the average LE M(t), where the av- 
eraging is performed over an ensemble of initial states 
or system Hamiltonians: M(t) ~ e~ At . The decay rate 
A equals the average Lyapunov exponent of the classi- 
cal system. This decay regime, known as the Lyapunov 
regime, provides a strong, appealing connection between 
classical and quantum chaos, and is supported by exten- 
sive numerical simulations J4|. For discussion of other 
decay regimes consult Ref. |2|. 

In this paper we report a new regime for the time de- 
cay of the unaveraged, individual LE for a semiclassical 
wave packet evolving in a two-dimensional billiard that 
is chaotic in the classical limit. We consider the general 
class of strong perturbations of the Hamiltonian that lo- 
cally modify the billiard's boundary: the perturbation 
only affects a boundary segment of length w small com- 
pared to the perimeter P, see Figs. Q] and [21 Both w and 
the perturbation length scale in the direction perpendic- 



ular to the boundary are considered to be much larger 
than the de Broglie wavelength A, so that the perturba- 
tion significantly modifies trajectories of the underlying 
classical system, see Fig. [TJ Our analytical calculations, 
confirmed by results of numerical simulations, show that 
the LE in such a system follows the exponential decay 
M(t) ~ e -27 ', with 7 being the rate at which classi- 
cal particles would escape from an open billiard obtained 
from the original, unperturbed billiard by removing the 
perturbation-affected boundary segment. The LE decay 
is independent of the shape of a particular boundary per- 
turbation, and only depends on the length of the per- 
turbation region. Furthermore, our numerical analysis 
shows that for certain choices of system parameters the 
exponential decay persists for times t even longer than 
the Heisenberg time in. 




FIG. 1: An unperturbed, chaotic billiard (solid-line), together 
with the perturbation (dashed-line) . The boundary of the 
unperturbed billiard consists of two segments, Bo and B\ . The 
perturbation replaces the latter segment by Bi , rendering the 
perturbed billiard to be bounded by Bo and B\. The initial 
Gaussian wave packet is centered at ro. Three possible types 
of trajectories, so, si and si, leading from ro to another point 
r, are shown. 



We proceed by considering a Gaussian wave packet, 



0o W = 



1 



■ exp 



i , s (r-r ) 2 



(2) 



centered at a point ro inside the domain A of a two- 
dimensional chaotic billiard (e.g. the solid-line boundary 
in Fig. [I}, and characterized by an average momentum po 
that defines the de Broglie wavelength of the moving par- 
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tide, A = ft/ 1 Po|- The dispersion a is assumed to be suffi- 
ciently small for the normalization integral J A dr\<fio(r)\ 2 
to be close to unity. We let the wave packet evolve in- 
side the billiard through a time t according to the time- 
dependent Schrodinger equation with hard-wall (Dirich- 
let) boundary conditions. This evolution yields the wave 
function 4>t{r) = {r\e~ lHt t h \4>o) , where H stands for the 
Hamiltonian of the billiard. Then, we consider a per- 
turbed billiard obtained from the original one by modi- 
fying the shape of a small segment of its boundary. Fig. Q] 
illustrates the perturbation: the unperturbed billiard is 
bounded by segments Bo and Bi, whereas the boundary 
of the perturbed billiard is composed of Bq and B\ . The 
perturbation, B\ — > B\, is assumed to be such that the 
domain A of the perturbed billiard entirely contains the 
domain A of the unperturbed one. The time evolution 
of the initial wave packet, Eq. ([2]), inside the perturbed 
billiard results to <f>t(r) = {r\e~ lHt / h \4>o) , with H being 
the Hamiltonian of the perturbed billiard. Then, the LE, 
defined in Eq. (JXJ) , reads 



M(t) = 



dr#(r)<Mr) 



(3) 



where the asterisk denotes complex conjugation. 

We now present a semiclassical calculation of the over- 
lap integral Eq. As the starting point we take the 
expression Q for the time evolution of the small (such 
that a is much smaller than the characteristic length scale 
of the billiard) Gaussian wave packet, defined by Eq. ((2]): 

0t(r) « (47R7 2 ) 1 / 2 Mr,r ,i)e-^-P°> 2 / 2n2 . 

s(r,r , t) 

(4) 

This expression is obtained by applying the semiclassi- 
cal Van Vleck propagator with the action linearized 
in the vicinity of the wave packet center ro, to the wave 
packet <po(r). Here, the sum goes over all possible tra- 
jectories s(r, ro, t) of a classical particle inside the unper- 
turbed billiard leading from the point ro to the point r 
in time t (e.g. trajectories sq and si in Fig. [Ij, and 



£«(r,r ,t) 



2irih 



exp 



1 a i ±\ • 7r ^ 
-S s (r,r ,i) - i — 



(5) 



where S s (r, ro,i) denotes the classical action along 
the path s. In a hard- wall billiard S s (r, ro,£) = 
(m/2£)Lg(r, ro), where L s (r, ro) is the length of the tra- 
jectory s, and m is the mass of the moving particle. In 
Eq. ®,.D, = | det(-0 2 S s /<9r<9r o )| is the Van Vleck de- 
terminant, and v s is an index equal to twice the number 
of collisions with the hard-wall billiard boundary that a 
particle, traveling along s, experiences during time t @. 
In Eq. (j4|), p s — —dS s (r,r ,t)/dr stands for the ini- 
tial momentum of a particle on the trajectory s. The 
expression for the time-dependent wave function 0*(r) 
is obtained from Eq. @ by replacing the trajectories 



s(r, r ,t) by paths s(r, ro,t), that lead from r to r in 
time t within the boundaries of the perturbed billiard 
(e.g. trajectories sq and S\ in Fig. [Ij. 

The wave functions of the unperturbed and perturbed 
billiards at a point r € A can be written as 



t (r)=«^(r)+^(r), 



(6) 



where <p[°\r) is given by Eq. ((4]) with the sum in the 
RHS involving only trajectories sq, which scatter only 
off the part of the boundary, Bq, that stays unaffected 
by the perturbation, see Fig. [TJ On the other hand, the 
wave function ( </>j ) involves only such trajectories 
si ( Si ) that undergo at least one collision with the 
perturbation-affected region, B\ ( B\ ), see Fig. [T] The 
LE integral in Eq. (|3|) has now four contributions: 



dr4>* t 4> t = 



[ dr 




2 + [ dr 
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[ dr 






C dr\ 
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id) 



(7) 



We argue that the dominant contribution to the LE over- 
lap comes from the first integral in the RHS of the last 
equation. Indeed, all the integrands in Eq. (|7|) contain 
the factor exp[i(S s — S s >)/H — i^{v s — ^«')/2], where the 
trajectory s is either of the type sq or si , and s' is either of 
the type sq or si, see Fig. [TJ An integral vanishes if there 
is no correlation between s and s', since the correspond- 
ing integrand is a rapidly oscillating function of r. This 
is indeed the case for the last two integrals: they involve 
such trajectory pairs (s, s') that s is of the type sq or si, 
and s' is of the type Si, so that the absence of correla- 
tions within such pairs is guaranteed by the fact that the 
scale of the boundary deformation is much larger than 
X. Then, we restrict ourselves to the diagonal approxi- 
mation, in which only the trajectory pairs with s = s 
survive the integration over r. The second integral in the 
RHS of Eq. ([7|) only contains the trajectory pairs of the 
type (so, si), and, therefore, vanishes in the diagonal ap- 
proximation. Thus, the only non-vanishing contribution 
reads 



f dr 
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<-/ 


'A 






J A 



dp 



-Pt(A) 



so 
TTh 2 



"^2 (P*o "Po) 



exp 



-£2 (P-Po) 



(8) 



where D SQ = | dct(9p So /9r)|, with p SQ being the initial 
momentum on the trajectory So(r, ro, t), serves as the Ja- 
cobian of the transformation from the space of final po- 
sitions r 6 A to the space of initial momenta p <E Vt (A) . 
Here, Vt(A) is the set of all momenta p such that a trajec- 
tory, starting from the phase-space point (ro,p), arrives 
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at a coordinate point r S A after the time t, while un- 
dergoing collisions only with the boundary Bq (and, thus, 
avoiding 2?i), see Fig. [U Thus, dr|c^ 0) | 2 is merely the 
probability that a classical particle, with the initial mo- 
mentum sampled from the Gaussian distribution, experi- 
ences no collisions with B\ during the time t. Therefore, 
if the boundary segment B\ is removed, this integral cor- 
responds to the survival probability of the classical par- 
ticle in the resulting open billiard. In chaotic billiards 
the survival probability decays exponentially [7| as e -7 *, 
with the escape rate 7 given by 



7 = i> 



w 



(9) 



where v — |po|/m is the particle's velocity, and A stands 
for the area of the billiard. Equation © assumes that 
the characteristic escape time 1 /j is much longer that the 
average free flight time t[. In chaotic billiards the latter 
is given by § tj = nA/vP, where P is the perimeter of 
the billiard. Condition tf <C 1/7 is equivalent to w -C P. 
In accordance with Eqs. ^ and the LE decays as 



M(t) ~ exp(-27t). 



Equation (| 10[) constitutes the central result of the paper. 
Together with Eq. ([9]) it shows that for a given billiard 
the LE merely depends on the length w of the bound- 
ary segment affected by the perturbation and on the de 
Broglie wavelength A = h/mv. It is independent of the 
shape and area of the boundary perturbation, as well as of 
the position, size and momentum direction of the initial 
wave packet. (We exclude initial conditions for which 
the wave packet interacts with the perturbation before 
having considerably explored the allowed phase space.) 

The decay rate 7, and thus the LE, are also related to 
classical properties of the chaotic set of periodic trajec- 
tories unaffected by the boundary perturbation, i.e. to 
properties of the chaotic repellor of the open billiard @ : 



7 = A r - h 



KS, 



(11) 



where A r is the average Lyapunov exponent of the re- 
pellor, and /iks is its Kolmogorov-Sinai entropy. Thus, 
Eqs. ((TDJ) and (fTTj) provide an interesting link between 
classical and quantum chaos. 

In order to verify the analytical predictions we sim- 
ulated the dynamics of a Gaussian wave packet inside 
a desymmetrized diamond (DD) billiard, defined as the 
fundamental domain of the area confined by four inter- 
secting disks centered at the vertices of a square. Accord- 
ing to the theorem of Ref. [l(| the DD billiard is chaotic 
in the classical limit. It can be characterized by the disk 
radius R, and the length L of the longest straight segment 
of the boundary, see Fig. [21 We consider the Hamiltonian 
perturbation that replaces a straight segment of length 
w of the boundary of the unperturbed billiard by an arc 
of radius r, see Fig. [2J In general, w < 2r. 




FIG. 2: Forward-time wave packet evolution in the unper- 
turbed DD billiard, followed by the reversed-time evolution 
in the perturbed billiard. The initial Gaussian wave packet is 
characterized by the size a = 12 and de Broglie wavelength 
A = 15/n; the arrow shows the momentum direction of the ini- 
tial wave packet. The DD billiard is characterized by L = 400 
and R = 200-v/lO. The perturbation is defined by w = 60 and 
r = 30. The propagation time corresponds to approximately 
10 collisions of the classical particle. 



To simulate the time evolution of the wave packet 
inside the billiard we utilize the Trotter- Suzuki algo- 



(10) rithm 11]. Figure [2] illustrates the time evolution of 



a Gaussian wave packet in the DD billiard followed by 
the time-reversed evolution inside the perturbed billiard. 
The parameters characterizing the system are L = 400, 
R = 200VT0, w = 60 and r = 30. The Gaussian wave 
packet is parametrized by a = 12, A = 15/7T; the ar- 
row shows the momentum direction of the initial wave 
packet. The evolution time t in Fig. [5] corresponds to 
some 10 free flight times of the corresponding classical 
particle, i.e. t = 10if . 

Figure [3] shows the time dependence of the LE com- 
puted for the DD billiard system characterized by L = 
400, R = 200\/T0, <7 = 3, and A = 5/tt. The initial mo- 
mentum direction is the same as in Fig. [2l The different 
LE decay curves correspond to different shapes of the lo- 
cal boundary perturbation: the width of the perturbation 
region stays fixed, w = 60, and the curvature radius of 
the perturbation arc takes the values r = 30; 35; 40; 45. 
In all four cases the LE displays the exponential decay 
for times t up to 40if-45if followed by LE fluctuations 
around a saturation value, M s . The thick solid straight 
line shows the trend of the e -27 * exponential decay, with 
7 given by Eq. ([9]). One can see strong agreement be- 
tween the numerical and analytical LE decay rates. We 
have also verified numerically that the LE decay rate 
is independent of the momentum direction of the initial 
wave packet. 

The inset in Fig. [3] presents the time decay of the av- 
erage LE M(t), with the averaging performed over 16 
individual decay curves M(t) corresponding to different 
values of the arc radius r, ranging from r = 30 to 45. 
The saturation mechanism for the LE decay was first 
proposed by Peres [l| and later discussed in Ref. Q . The 
LE saturates at a value M s inversely proportional to the 
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FIG. 3: (Color online) The Loschmidt echo (LE) decay in the 
DD billiard for four different values of the curvature radius 
r of the arc perturbation. The width of the perturbation 
region is fixed, w = 60. The other system parameters are 
L = 400, R = 200 VW, a = 3, and A = 5/tt. The solid straight 
line gives the trend of the exp(— 2^t) decay, with 7 given by 
Eq. ((9)1 . The inset presents the decay of the average LE, with 
averaging performed over individual LE curves corresponding 
to different values of r. 



number N of energy levels significantly represented in 
the initial state. If the areas of the unperturbed and 
perturbed billiards are relatively close, then N sa A /a 2 
and M s ~ a 2 /A. (We have verified the latter relation 
by computing the LE saturation value for billiards of dif- 
ferent area.) Thus, one might expect the exponential 
decay of the LE to persist for times f < t s , with the 
saturation time t s = (1/271;) In iV. The latter can be 
longer than the Heisenberg time £h = A/2-kXv for a sys- 
tem with sufficiently large effective Hilbert space, since 
ts/tn ~ (X/w)\nN. Indeed, for the system correspond- 
ing to Fig. [3] one has in ~ 29if, whereas the exponential 
decay persists for times t < 40if. 

Finally, we sketch a principal experimental scheme 
for measuring the LE decay regime proposed in this 
paper. Consider a two-dimensional, AlGaAs-GaAs 
heterojunction-based ballistic cavity with the shape of a 
chaotic billiard, e.g. Fig. [TJ Let the initial electron state 
to be given by \^o) = |0o) ® |x), where \<po) is the spatial 
part defined by Eq. J2j , and \ X ) = (||) + |1)) rep- 

resents a spin-1/2 state. Here, 1 1) and 1 1) are the eigen- 
states of the spin-projection operator in the z-direction, 
perpendicular to the billiard plane. Then, \x) is the 
eigenstate of the spin-projection operator s x = ^cr x , with 
a x = I T)(4 I + I l)(T |j m some cc-direction, fixed in the 
billiard plane. Suppose now that a half-metallic ferro- 
magnet, magnetized in the z-direction, is attached to the 
boundary of the ballistic cavity. (One may consider the 
region bounded by Bq and B\ in Fig. Q] to represent the 
ballistic cavity, and the region bounded by B\ and B\ 
to represent the ferromagnet.) Then the ferromagnet- 



cavity interface will reflect the | t)-component of the 
state, but will transmit the | |)-component. As a re- 
sult, the two components will evolve under two different 
spatial Hamiltonians, H and H, corresponding to the ge- 
ometry of the ballistic cavity and the geometry of the 
cavity- ferromagnet compound respectively. Then |\I>o) 
will evolve to 

|* t ) - -J= \e-' lHt ' h \^) ® I T) + e-^/^o) ® 1 1) 

(12) 

The expectation value of the projection of the spin in the 
x-direction is related to the LE overlap by 



s (t) = (* t |«x|*t) 



-Re( 



b \e im ' h e~' lHt l h 



\<M, (13) 



where Re denotes the real part. As we have shown 
above, this overlap is real and decays exponentially in 
time. Therefore, the average spin projection in the re- 
direction will also relax exponentially with time, i.e. 
s x (t) ~ 4exp(— •yt), with the relaxation rate 7 deter- 
mined by Eq. @. This result provides a link between 
the spin relaxation in chaotic, mesoscopic structures [l2| 
and the LE decay due to local boundary perturbations. 
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